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ABSTRACT 

Jj 

The  generalized  likelihood  ratio  detector,  for  deciding  between  the  right- 

left  ambiguity  of  a  line  array  attempting  to  estimate  the  angle  of  arrival  of 

a  plane  wave,  is  derived.  Two  scenarios  are  considered,  the  first  with  noisy 

measured  antenna  angle,  the  second  with  noiseless  antenna  angle  measurements. 

The  detector  for  both  cases  is  a  cross-correlator  of  the  sample  ac  components 

of  the  measured  antenna  and  source  angle  waveforms.  ^ 

r 
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INTRODUCTION 


A  line  array  inherently  has  a  cone  of  ambiguity  in  its  response.  When 
the  array  lies  in  the  horizontal  plane,  and  a  source  is  located  in  that  same 
plane,  the  ambiguity  reduces  to  a  right-left  uncertainty,  which  cannot  be 
resolved  without  some  manuevering  on  the  part  of  the  source  or  array.  If  the 
line  array  is  moving  randomly,  unintentionally  or  uncontrollably,  this 
movement  can  serve  as  a  means  of  making  a  high  quality  decision  about  the 
source  direction,  if  the  array  angle,  as  well  as  the  source  angle  relative  to 
the  line  array,  are  measured. 
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Figure  1.  Geometry  of  Line  Array  and  Source 

The  situation  of  interest  here  is  described  in  figure  1.  The  line  array 
is  being  towed  due  north;  however,  it  is  undergoing  rigid  bar  rotation  about 
this  direction  in  a  random  manner,  as  described  by  random  process  a(t),  which 
is  the  actual  antenna  angle  relative  to  the  towing  direction. 

The  actual  source  angle,  relative  to  the  line  array  end-fire  direction, 
is  s(t).  Furthermore,  the  actual  source  angle,  relative  to  the  towing 
direction,  is  e,  an  unknown  constant;  it  is  presumed  that  »  is  constant 
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throughout  the  observation  interval.  Reference  to  figure  1  reveals  that  these 
various  quantities  are  interrelated  according  to  the  alternatives 


p(t)  +  s(t)  for  hypothesis  1, 
|a(t)  -  s(t)  for  hypothesis  2,  H£ 


(1) 


However,  it  is  unknown  which  hypothesis  is  correct;  nevertheless,  it  is 
desired  to  make  a  reliable  decision,  so  that  an  accurate  estimate  of  the 
source  direction  can  be  made.  From  (1),  observe  that  we  can  express 


s(t) 


f"o  -  a(t)  for  Hj^ 
t-e  +  a(t)  for  H2 


(2) 


which  will  be  needed  in  later  developments. 
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NOISY  MEASURED  ANTENNA  ANGLE 


In  this  section,  the  measured  antenna  angle  is  not  a(t)  as  desired,  but 
rather  is 

x ( t )  *  a(t)  +  m(t)  ,  (3) 

where  m(t)  is  an  unavoidable  additive  noise  process.  Also  the  measured  source 
angle  is  not  s(t),  but  instead  is 

y(t)  «  s(t)  +  n(t)  ,  (4) 

where  n(t)  is  likewise  an  undesirable  additive  perturbation,  due  to  limited 
observation  time,  array  length,  ambient  noise,  etc.  The  three  random 
processes  a(t),  m(t),  n(t)  are  presumed  to  be  zero-mean  Gaussian  processes, 
independent  of  each  other. 

Combining  (2)-(4),  the  situation  is  as  follows:  the  available 
measurements  upon  which  a  decision  must  be  reached  are  the  two  waveforms 


x { t )  =  a(t)  +  m(t) 


y(t)  -  o  -  a(t)  +  n(t) 


for  Hi 


x(t)  =  a(t)  +  m(t) 


y(t)  =  -9  +  a(t)  +  n(t) 


for  H, 


On  the  basis  of  waveforms  x(t)  and  y(t),  what  is  the  best  decision  and  what  is 
the  corresponding  estimate  of  9? 
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Derivation  of  Generalized  Likelihood  Ratio 


Let  a  be  the  time  sampling  increment  applied  to  measurement  waveforms 
x(t)  and  y(t);  assume  that  the  samples  of  the  three  processes  are 
statistically  independent  at  this  rate.  Denote 

\  *  x(ka)  for  1  <  k  <  K  , 

(6) 

y^  *  y(ka)  for  1  <  k  <  K  , 

wnere  Ka  is  the  total  observation  time,  and  let  the  collections  of  samples  be 
denoted  by 

X  =  xv  x2,  ....  x^ ,  Y  «  yj,  y2,  ....  yK,  A  =.  ax,  a2 . aK  .  (7) 

Tnen  for  a  fixed  A  and  a  hypothesized  value  ©h  for  the  source  angle,  the 
conditional  probability  density  function  under  H^,  of  the  total  set  of 
measurements  X,  Y,  is 


pl(x,  r|v  .  i{L  L  k-^j)  . 


wnere  the  normalized  Gaussian  probability  density  function  is 

*(t)  -  (2w)‘i/2  exp(-t2/2)  .  (9) 

Here  we  used  (5A)  and  the  Gaussian  character  of  processes  m(t)  and  n(t)  with 
standard  deviations  om  and  on  respectively. 

We  now  must  weight  (8)  by  the  Gaussian  probability  density  function  for 
process  a(t)  and  integrate  over  A,  to  determine  the  unconditional  probability 
density  function  of  X,  Y,  for  hypothesized  value  0^  Carrying  out  the 
integrals  and  simplifying  the  result,  we  obtain 
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Pit*.  yK>  -  (2'Wn  *\'2) 


-K 


exp 


ERan  r  2 


am  r,  „  .2 

i^|(yk  11  *  r3. 


1 

Vn  "  J 


where  we  define 


'am 


1  +  1  „  1.1  „  1+1.1 
“7  +  “7  ’  Ran  *  "T  +  “T  ’  R3  =  77  +  77  +  “T  > 


0.  o* 
a  m 


o“  o„ 
a  n 


0“  o-  o*- 

a  m  n 


(10) 


(11) 


and  use  the  shorthand  notation 


_  K 

2-2  •  <12> 

k  k-1 

In  a  similar  fashion,  the  unconditional  probability  density  function  of  X,  Y 
under  H2,  for  hypothesized  value  ©h,  is  given  by 


P2«.  *K>  ■  (z-.Vn  «3/2) 


-K 


*  exp 


an 


2R3°m  H 


z  24- 


am 


2R3a 


7  2(yk+®h^  +  — 1  7-7  HMVVr 
n  ^  R3aman  *  J 


(13) 


Now  if  ^  were  known,  we  could  evaluate  the  likelihood  ratio  by  taking 
the  ratio  of  (10)  and  (13).  However,  we  must  resort  instead  to  a  generalized 
likelihood  ratio,  by  computing  the  two  values  of  ©h  that  maximize  (10)  and 
(13)  respectively,  and  then  taking  the  ratio  of  the  two  maxima  [1,  p.  92]. 
This  procedure  is  not  optimum  in  any  sense;  however,  it  often  leads  to 
physical  processors  that  perform  well. 


The  values  of  ©u  that  maximize  (10)  and  (13)  are  given  respectively  by 
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a2 

*2  "  +  "T“T  K  fxk 

K  a,  +  <J_  n 

a  m 

These  results  have  a  reasonable  physical  interpretation:  From  (5A),  the  sum 
x(t)  +  y(t)  would  eliminate  the  random  process  a(t),  and  the  sample  mean  of 
the  sum  would  give  an  estimate  of  e  under  Hj.  However,  m(t)  contaminates 

2  2  2 

the  a(t)  contribution  according  to  (3);  thus  the  scale  factor  ^/(a*  +  o^)  in 

(14A)  indicates  how  trustworthy  the  sample  mean  of  x(t)  is.  The  noise  n(t)  in 
y(t)  is  unavoidable  but  is  partially  suppressed  by  the  inherent  averaging  of 
the  sample  mean  of  y(t).  A  similar  argument  holds  for  x(t)  -  y(t)  under  H2. 

The  logarithm  of  the  generalized  likelihood  ratio  is  (proportional  to) 


for  Hj  ,  (14A) 


for  H2  ,  (14B) 


+ 


1 

7 


^am  °m 


+  7^xk^k  “  »l>  * 


(15) 


where  we  used  (10)  and  (13)  with  and  ©2  substituted  for  Now  let 

sx-|xk-  Sy  ■  '  P-£Vk-  (: 

Then  (15)  becomes,  upon  use  of  (14)  and  simplification  and  cancellation  of 
various  terms, 


l>.  It, 


■  l»> 
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T  -  P  -  K 


■£**  -  rfx»2n 


(17A) 


where 


(17B) 


xk  “  xk  ~ 

Wi 


■t5>,  forl<k<K, 


»k-yk-Z?,n  for  1  i  k  i  K  ■ 


are  defined  as  the  sample  ac  components  of  measurements  X  and  Y;  that  is,  the 
sample  means  are  subtracted  from  the  measurements. 

The  generalized  likelihood  ratio  test  in  (17)  says  to  cross-correlate  the 
sample  ac  components  of  both  measured  waveforms  and  to  compare  with  zero 
(assuming  Hj  and  H2  are  equally  likely  apriori).  That  is,  the  test  is 


K _ H? 

T  =  xkyk  $  0  ‘ 
k=l  H, 


Observe  that  this  decision  rule  makes  no  use  of  the  variances  of  any  of  the 
processes  a(t),  m(t),  n(t),  although  this  information  was  presumed  known  in 
the  above  derivation.  Of  course,  the  source  angle  estimates  in  (14)  do 

2  2 

require  knowledge  of  the  signal-to-noise  ratio  oa/om  in  the  x(t)  measurement 

of  the  antenna  angle  a(t).  Once  the  decision  of  Hj  yS  H2  is  made  via  (19), 
the  corresponding  estimate  of  the  actual  source  angle  ©  is  taken  from  (14). 
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On  the  Performance  of  Test  (19) 


From  (5)  and  (6),  we  find  that 


Vk 


for 

►  * 

for  Hg 


xmyn  “  0  for  m  *  n  •  ( 20 ) 

Then  (17A)  yields  the  mean  value  of  the  generalized  likelihood  ratio  test 
statistic  as 


r-u-i)  o2a 

I  4 


(21) 


Since  the  statistics  of  T  are  desired  different  under  the  two  hypotheses, 

(21)  indicates  that  large  K  and  are  desired.  That  is,  a  large 
observation  time  and  a  widely-moving  antenna  give  better  performance  of  the 
test;  both  of  these  conclusions  are  physically  plausible. 


Although  generalized  likelihood  ratio  test  (19)  does  not  require 
knowledge  of  any  variances,  the  performance  (in  terms  of  the  error 
probability)  does  depend  on  all  the  variances.  However,  the  performance  does 
not  depend  on  the  actual  value  e  of  the  source  angle.  To  see  this,  we  employ 
(5)  in  (18)  to  obtain 


-ak  *  nk  -  (-Sa  +  Sn)/K 


ak  +  nk  -  <  Sa  +  Sn>'K 


for  H 


for  H„ 


i 


(22) 
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sa  - 2\-  s„  ’1\  ■  <23> 

k  x 

Thus  o  is  absent  from  (22),  and  since  o  is  not  involved  in  x(t)  or  xk,  test 
statistic  T  is  independent  of  e. 

If  we  develop  (19)  in  more  detail  and  make  use  of  (18)  and  (5),  we  find 
we  can  express 


Ti  -  *  mk><\  -  nk>  -  k  £<ak  *  "k»ak  -  nkjJ 

T2  ■  [2<ak  *  mk><ak  *  "k>  -  1  ?<ak  *  mk>(ak  *  nk] 


for  , 


for  H2  . 


Tnus  the  statistics  of  are  identical  to  those  of  -T2. 


(24) 


Based  upon  the  results  in  [2,3],  an  exact  analysis  of  the  cumulative  and 
exceedance  distribution  functions  of  test  statistic  (17A)  is  possible  and  will 
be  documented  in  a  NUSC  technical  report  shortly;  in  fact,  a  more  general 
processor,  where  the  sample  means  term  is  scaled  prior  to  subtraction,  will  be 
analyzed. 
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NOISELESS  MEASURED  ANTENNA  ANGLE 

In  this  section,  the  noise  m(t)  in  the  antenna  angle  measurement  is  zero; 
thus,  from  (3), 


x(t)  =  a(t)  (25) 

under  and  H2.  We  also  remove  the  Gaussian  assumption  on  the  statistics 
of  antenna  movement  a(t),  and  allow  a(t)  to  be  completely  general. 
Furthermore,  we  allow  any  statistical  dependence  among  the  samples  A  of  a(t) 
in  (7).  However,  we  retain  the  Gaussian  assumption  on  the  additive  noise 
process  n(t)  in  (4)  and  (5),  and  keep  the  statistical  independence  of  its 
samples  (n(ka)} 

Derivation  of  Generalized  Likelihood  Ratio 


Based  upon  these  premises,  the  conditional  probability  density  function 
under  Hj  0f  measurements  X,  Y,  for  a  fixed  A  and  hypothesized  eh,  is 

Pl(X,  Ylv  A)  =  jr  f.K  -  ak)  I-  *  ■-*)]  •  (26) 

(This  is  also  the  limit  of  (8)  as  om-*o+.)  Then  letting  joint  probability 
density  function  pa(A)  represent  the  arbitrary  statistical  dependence  of 
samples  A,  the  unconditional  probability  density  function  under  of  X,  Y  is 

Pl(X,  Y J©h)  =  JdA  pa(A)  Pl(X,  Y |®h.  A)  = 

using  (26),  (9),  and  the  sifting  property  of  delta  functions.  (If  measurement 
noise  m(t)  in  (3)  were  non-zero,  this  simplification  of  the  probability 
density  function  in  (27)  would  not  be  possible.) 
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The  value  of  ^  that  maximizes  probability  density  function  (27), 
regardless  of  the  form  of  the  probability  density  function  pa>  is 

•l  *  7^{xk  +  yk>  *  (28) 

which  is  simply  the  sample  mean  of  the  sum  waveform  x(t)  +  y(t).  This  result 
is  consistent  with  the  earlier  one  in  (14A),  for  am  .  o,  and  the  ensuing 
discussion. 

In  a  similar  fashion,  the  probability  density  function  of  X,  Y  under  H2 
is 

P2(X,  V|%)  .  <f7T'n)-*  expjl  Jy  %iyk  *  ph  -  xk)2]  p,(X)  ,  (29) 

and  the  maximizing  choice  of  ©h  is 

®2  *  ^  ~  ^k^  *  (30) 


which  is  the  sample  mean  of  difference  waveform  x(t)  -  y(t). 


There  follows,  from  (27)-(30),  the  logarithm  of  the  generalized 
likelihood  ratio  as 


2/i  y 

2ffn pxix,  yjoj) 


-l(yk  *  ®2  -  \)z  +2(yk 


(31) 


The  generalized  likelihood  ratio  test  for  the  two  hypotheses  in  (1)  is 
therefore 
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or  up on  substitution  of  (28)  and  (30),  and  use  of  (18),  simply 


Sly.  >2  0  . 

k  k  < 
k  Hi 


x(t)  +  y(t)  »  »  +  n(t)  under  Hj  , 

x{ t)  -  y(t)  *  o  -  n(t)  under  H2  .  ( 

That  is,  except  for  zero-mean  measurement  noise  n(t),  the  sum  waveform  is 
constant  under  Hj,  whereas  the  difference  waveform  is  constant  under  H£. 

An  alternative  form  for  (32)  and  (33)  is 


msA 


?<xk  +  *k  '  *1>2  I2  Sxk  -  *k  -  *2)Z 
*  Hx  * 


As  in  tne  previous  section,  the  cross-correlation  of  the  sample  ac  components 
of  the  measurements  X,  Y  should  be  compared  with  zero.  This  decision  rule 
holds  for  any  statistics  of  antenna  movement  a(t). 

The  alternative  form  in  (32)  has  an  interesting  interpretation: 

Reference  to  (27)— (28)  reveals  that  ©j  is  the  best  constant  fit  to 
{xk  +  y|<}5  in  a  least  squares  sense.  Thus  the  left  side  of  (32)  is 
the  actual  value  of  the  least  squares  error  of  a  constant  fit  to  the  sum 
waveform.  Similarly,  the  right  side  of  (32)  is  the  least  squares  error  of  a 
constant  fit  to  the  difference  waveform.  Whichever  error  is  smaller,  that 
hypothesis  is  selected.  This  decision  rule,  (32),  is  consistent  with  the 
observation  from  (5)  and  (25)  that 


m 


'■V«\YVAV.NV 
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T  *  ?Vk  -  K  |xm  2*n  l2  «. 


*  * 


just  as  in  (17A) .  All  the  ensuing  discussion  there  through  (24)  is  directly 
relevant  for  this  case  as  well. 
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SUMMARY 


The  generalized  likelihood  ratio  test  statistic,  for  both  the  noisy  as 
well  as  the  noiseless  antenna  angle  measurement,  is  a  cross-correlator  of  the 
sample  ac  components  of  the  measured  antenna  and  source  angle  waveforms. 

Exact  performance  of  this  processor  can  be  accomplished,  since  the  test 
statistic  is  a  quadratic  form  of  correlated  Gaussian  random  variables;  in 
fact,  the  complete  cumulative  and  exceedance  distribution  functions  of  the 
quantity 


K  K  K 

2  V*  -  2  yn 

k=»l  rrwl  Hal 


(36) 


for  any  scaling  y,  is  capable  of  exact  analysis  and  will  be  presented  in  a 
future  NUSC  technical  report. 
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